HST /COS SPECTRA OF DF TAU AND V4046 SGR: FIRST 
DETECTION OF MOLECULAR HYDROGEN ABSORPTION 
AGAINST THE Lya EMISSION LINE ^ 



Hao Yang 

JILA, University of Colorado and NIST, Boulder, CO 80309-0440 
haoyangSj ilaul . Colorado . edu 
Jeffrey L. Linsky 
JILA, University of Colorado and NIST, Boulder, CO 80309-0440 
j linskyO j ilaul . Colorado . edu 
Kevin France 

CASA, University of Colorado, Boulder, CO 80309-0389 
kevin . f ranceOcolorado . edu 

ABSTRACT 

We report the first detection of molecular hydrogen (H2) absorption in the 
Lyman-a emission line profiles of two classical T Tauri stars (CTTSs), DF Tau 
and V4046 Sgr, observed by HST /COS. This absorption is the energy source for 
many of the Lyman-band H2 fluorescent lines commonly seen in the far-ultraviolet 
spectra of CTTSs. We find that the absorbed energy in the H2 pumping transi- 
tions from a portion of the Lyman-o; line significantly differ from the amount of 
energy in the resulting fluorescent emission. By assuming additional absorption 
in the H I Lyman-a profile along our light of sight, we can correct the H2 absorp- 
tion/emission ratios so that they are close to unity. The required H I absorption 
for DF Tau is at a velocity close to the radial velocity of the star, consistent with 
H I absorption in the edge-on disk and interstellar medium. For V4046 Sgr, a 
nearly face-on system, the required absorption is between +100 km s~^ and +290 
km s~^, most likely resulting from H I gas in the accretion columns falling onto 
the star. 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Science Institute. STScI is operated by the Association of Universities for 
Research in Astronomy, Inc. under NASA contract NAS 5-26555. 
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Subject headings: stars: pre-main sequence — ultraviolet: stars — accretion, 
accretion disks — stars: individual (DF Tau, V4046 Sgr) 



INTRODUCTION 



Molecular hydrogen (H2) emission lines are commonly observed in the far-ultraviolet 
(FUy ) spectra of classical T Tauri stars (CTTSs) Jflerczeg et al.ll2002l . [20oi . l20od : lArdila et al, 
2OO2I ) ■ iBrown et al.l (Il98ll ) first detected these lines fr om T Tau using lU E. As suggested by 
detections of identical lines in the sunspot spectrum (jjordan et al.l 119771 ) , these fluorescent 
H2 lines are thought to be photoexcited from the ground electronic state to the B (or C) 
electronic state primarily by coincidence with hydrogen Lyman-a but also by other strong 
atomic emission lines (e.g., C II, C IV, and O VI ) in the UV. 

The H2 fluorescence may arise from various locations in protostellar sys tems. In their 



analys is of HST STIS E140M spectra of 6 pre-main sequence (PMS) stars, iHerczeg et al. 



(I2OO6I ) found that blueshifted H2 lines of RU Lupi, T Tau, and DG Tau are likely formed in 
outflows, while the H2 lines of TW Hya, DF Tau, and V836 Tau show no radial velocity shifts 
from the photospheric lines and are likely formed in warm (~2500 K) surfaces of their cir- 
cumstellar disks. Stars such as T Tau also show spatially extended H2 fluorescent emission in 
associated nebulosit y, which is likely pu mped by local shocks and outflows rather than stellar 
Lyman-a emission (IWalter et al.ll2003l ). For the diskless counterparts of C TTSs, the naked 



T Ta uri stars (NTTSs), H2 features are not seen in their FUV spectra (jingleby. L. et al 



2009 



Yang et al. 2011, in preparation), indicating that H2 fluorescent emission requires the 
presence of H2 gas close to the central star. Studying the H2 fluorescent emission therefore 
provides valuable information on the physical properties of protoplanetary disks, which are 
90% composed of H2 gas. 

While the H2 fluorescent lines have been studied in a number of CTTSs, the pumping 
transitions had not been observed in absorption against the Lyman-a emission line. In 
this Letter, we present new FUV spectroscopy of two CTTSs, DF Tau and V4046 Sgr, for 
which we detect such absorption for the flrst time as a result of the very low noise and high 
throughput of the Cosmic Origins Spectrograph (COS) on HST. 



2. OBSERVATIONS AND DATA REDUCTION 



We observed DF Tau (RA = 04:27:02.795, DEC = 57:12:35.38) and y4046 Sgr (RA 
= 18:14:10.466, DEC = -32:47:34.50) with the COS JdIxou et al.l I2OI0I : lOsterman et"aL 
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20 111 ) on January 11 and April 27, 2010, respectively. COS is a high-throughput, moderate- 
resolution UV spectrograph installed on the HST in May 2009. During our HST GTO 
program 11533, we used both the G130M and G160M gratings of the COS FUV channel to 
cover the 1136 A-1796 A region. Since there is a small gap (~15 A in wavelength coverage) 
between the two segments of the COS detector, we observed each star at 4 central wavelength 
settings for both gratings to provide continuous spectral coverage and minimize any fixed- 
pattern noise. The total exposure time for each star was about 10,000 sec during 4 HST 
orbits. The spectral resolution was approximately 17, 000-18000, with extended wings in 
the line-spread function. The exten ded wings are induced by polishing errors on the HST 
primary and secondary mirrors (see iGhavamian et al.ll2010l ). 



We reduced the DF Tau and V4046 Sgr spectra using the COS calibration pipeline, 
CALCOsE(v 2.12, March 19, 2 010), and combined them with a custom IDL coaddition routine 
described by iDanforth et al.l (120101 ). In Figure [H we show two portions of the coadded FUV 
spectra for both stars as examples. To match the atomic and molecular emission lines with 
laboratory wavelengths, we corrected for the radial velocity of DF Tau (+15 km s~^), but 
found that an additional correction of -|-8 km s~^ is needed because of inaccuracies in the 
CALCOS wavelength solution. The radial velocity of V4046 Sgr is close to 0, and the 
wavelength solution is accurate. 



3. DF TAU AND V4046 SGR 



DF Tau is a binary system consisti ng of two early M stars separated bv 0.1^^ (jSchaefer et al 



2mm- Its disk is inclined by 80-85° flJohns-KruU fc Valentill200ll : lArdila et al.ll2002f l. i.e., 
the disk is seen nearly edge-on. The dist ance to DF Tau is ge nerally adopted as 140 pc, the 
distance to the Taurus Molecular Cloud. iHerczeg et al.l (120061 ) analy zed in detail the H2 flu- 



oresce nt lines of DF Tau observed with the STIS E140M grating. IHerczeg &: Hillenbrand 
(I2OO8I ) estimated a visual extinction {Ay) of 0.6 mag and accretion rate in the range of 
2.3-4.6 X 10-^Mq yr-i for DF Tau. 

V404 6 Sgr is a spectroscopic b inary with a separation of 9 Rq and an orbi tal period of 



2.42 days (IStempels fc Gahmll2004l ). The pair consists of a K7V and a K5V star ( lOuast et al. 
2OOOI ) . At a distance of ~ 72 pc (Tor res et al. 2008), it is a n isolated system and the 
extinction is pr actically Ay = 0.0 mag (IStempels fc Gahml |200J) . The circumbinary disk is 



inclined at 35° ( Ouast et al. 2000; Kastner et al. 20081 ). somewhat face-on. V4046 Sgr may 



*We refer the reader to the Cycle 18 COS Instrument Handbook for more details: 
http : //www . stsci . edu/hst/cos/documents/handbooks/current/ cos_cover . html) 
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be a member of the Pic Moving group and could be as old as 12 Myr (Ortega et al. 2002). 



Jensen &: Mathieul ( 119971 ) found that it has little excess emission in the near-infrared (near- 



IR) wavelengths but a large excess emission at longer wavelengths, suggesting that optically 
thick dust in the regions of the disk close to the star has been cleared out. We know of no 
previous detailed study of H2 emission from V4046 Sgr. 

Besides the H2 line emission studied in this work, the FUV continuum, which is also 
detected in DF Tau and V4046 Sgr, will be analyzed in detail by K. France et al. (2011, 
submitted). 



4. ANALYSIS & RESULTS 



The FUV spectra of DF Tau and V4046 Sgr are dominated by H2 fluorescent lines 
(see Figure [T]). Below 1200 A, there are also many Werner-band H2 lines, which will be 
described in a separate paper (Yang et al. in preparation). In Figure [21 we show the 
Lyman-a profiles of DF Tau and V4046 Sgr. The apparent absorption features seen against 
the Lyman-a emission line have depths much greater than the noise at those wavelengths 
and coincide in wav elength with the Lyman-band H2 pumping transitions (see Table 3 of 
Herczeg et al.ll2006l ). We first mask out the absorption features and fit 5th or 6th order 
polynomial curves to the Lyman-a profiles. To measure the amounts of absorbed energy, 
we integrate the area between the fitted curves and the observed spectra. The uncertainties 
in such measurements are mainly caused by the somewhat subjective determinations of the 
Lyman-a profiles without absorption. Some features result from the blended absorption of 
two or three pumping transitions, and the absorbed energy for these features represents the 
total absorption of the transitions. For DF Tau, 5 features are measured that correspond 
to 9 pumping transitions: 1218.52 A + 1218.57 A, 1219.09 A + 1219.10 A + 1219.15, 
1219.37 A + 1219.48 A, 1219.74 A, and 1220.18 A. For V4046 Sgr, 8 features are measured 
that correspond to 11 pumping transitions: 1212.43 A + 1212.54 A, 1213.36 A, 1213.68 A, 
1217.03 A + 1217.04 A, 1217.20 A, 1217.41 A, 1217.64 A, and 1217.90 A + 1217.98 A. 

We detect between 2 and 19 fluorescent H2 lines in the progressions produced by each 
pumping transition observed in absorptio n against the Lyman -o; lines. The H2 fluorescent 



lines are identified based on the line list of lAbgrall et al.l f 



we used a custom IDL fitting procedure (iFrance et al 



1993). To measure the line fluxes. 



2OIOI ) that convolves a Gaussian 



profile with the COS line spread function (LSF) to fit the observed H2 line profiles. The 
convolution of a Gaussian profile with the COS LSF only changes the shape of the profile 
but not the total line fiux. The uncertainties in the line fiuxes are generally less than 5 % 
for unblended lines, indicative of the high signal-to-noise of the data, typically > 50 for the 
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H2 lines. We next convert the line fluxes in each progression to the total energy emitted 
from the pumped upper level. Each H2 line in a given progression yields an estimate of the 
total energy in the upper level from the line theoretical branching ratios. We average the 
estimated total energy emitted from each upper level using only the strong unblended lines, 
and their standard deviations are < 15% of the mean values. Th e emission from each upper 



level is corrected for its dissociation probabihty as calculated by lAbgrall et al.l (I2OOOI ). The 



dissociation probability is typically zero or only a few percent. The absorptio n and emission 



fluxes for DF Tau are also corrected for extinction using Ay = 0.6 mag and the lCardelli et al. 
fll989h extinction law. 



For each absorption feature observed against the Lyman-a emission line, we have es- 
timated the absorbed energy (Eats) and the emitted energy (-Eem) from the corresponding 
H2 upper level. If there are no additional sources of absorption or emission, then the global 
average of Eats/Eem should be unity. We begin the analysis by assuming that Eabs/Eem along 
our line of sight should be close to unity, but the measurements show otherwise. We do flnd 
that the three features in the Lyman-a blue wing for V4046 Sgr have Eabs/Eem within a 
factor of two of unity (1.38, 0.99, 0.56), but the Eabs/Eem ratios at longer wavelengths are 
smaller than one by factors of 4-60 (see the fourth column of Tabled]). We therefore pro- 
pose that additional Lyman-a absorption in the line of sight between the location where 
H2 is pumped and the observer has reduced the observed absorbed energy in the pumping 
line and thus the Eabs/Eem ratios. To model the additional hydrogen absorption, we calcu- 
late Lyman-a absorption profiles using the Voigt function for a range of hydrogen column 
densities (Nhi) and radial velocities (vrad)- Then for each combination of Nhi and Vrad, we 
calculate the optical depths, t(A), at the wavelengths of the pumping lines. We correct the 
observed Eabs/Eem ratios by multiplying by e"'^^'*'^. For V4046 Sgr, the absorption features 
at 1217.03 A and 1217.20 A are not used. The 1217.03 A feature is close to the line center 
of Lyman-a, and the measurement is greatly affected by the absorption in the line center. 
The 1217.20 A progression does not have enough strong unblended emission lines to provide 
an accurate estimate of the emission energy. 

In Figure [3l we show for DF Tau and V4046 Sgr the combinations of Nhi and Vrad that 
can correct the Eabs/Eem ratios to be close to unity. For DF Tau, if we assume that the 
extinction is all due to interstellar dust, we can convert Ay = .6 mag to a total hydroge n 



column density of \og{NH) = 21.03 according to the relation in iPredehl fc SchmittI (119951 ). 
Note that this value assumes a standard gas-to-dust ratio for the interstellar medium. We 
mark in the top p anel of Figure the corresponding Nhi values assuming 100% and 50% 



neutral hydrogen. iHerczeg et al.l (120061 ) estimated for DF Tau the absorption against the 
red wing of the Lyman-a emission line and measured a log(NHi) of 20.75. This value is close 
to the hydrogen column density converted from Ay with 50% neutral content and represents 
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the neutral hydrogen in both the interstellar medium along the line of sight and possibly 
the edge-on disk of the system. As shown in Figure [3l within reasonable ranges of Nhi, the 
additional absorption required to bring all of the Eahs/Eem ratios close to unity for DF Tau 
requires radial velocities close to zero. Given that the uncertainty in Ay could be as large 
as 0.4 mag, and the ionization fraction and the gas-to-dust ratio in the disk are unknown, 
we think that the absorption is likely caused by a combination of interstellar medium and 
neutral hydrogen in the edge-on disk, though absorption by the H I columns in stellar winds 
or accretion columns can not be completely ruled out. 

For V4046 Sgr, which suffers negligible continuum extinction, the absorption required 
has a radial velocity between 100 and 290 km s~^, as shown in the bottom panel of Figure [31 



If Ay is close to zero (IStempels &: Gahnj|2004l ). then the radial velocity is close to 290 km s ^ . 



For this case, we list in Tabled] the calculated optical depth for the additional absorption as 
well as the measured and corrected Eats /E em for the absorption features detected in V4046 
Sgr. This is consistent with a scenario in which H I in the accretion columns is absorbing the 
red wings of the Lyman-a emission line in this system, which is oriented somewhat face-on. 



5. DISCUSSION 

Thanks to the excellent sensitivity and low background of COS, we were able to detect 
for the first time H2 absorption against the Lyman-a emission line profiles in two CTTSs. 
Because of the large aperture of COS (2.5"), the center of the Lyman-a line is filled with 
geocoronal Lyman- g emission and not u sable (see Figured]). The STIS E140M spectrum of 



DF Tau reported by lHerczeg et al.l (120061 ) shows that the line center and blue wing of Lyman- 
a are completely absorbed. Interstellar absorption must be responsible for the disappearance 
of the line center, because we see many H2 lines pumped by transitions coincident with the 
center of Lyman-a. On the other hand, we detect only a few H2 lines pumped by the 
transitions blueward of line center, suggesting that the blue-wing emission of Lyman-a has 
been absorbed by the stellar wind before the blue-wing radiation reaches the molecular gas 
in the disk. Since the disk of DF Tau is viewed nearly edge-on, the stellar wind must be 
present near the stellar equator to absorb the blue wing of the Lyman-a emission line. 

Understanding the geometry of the V4046 Sgr system requires more detailed considera- 
tion. The stellar wind must be weak for this somewhat face-on system since the blue wing of 
Lyman-a is not totally absorbed. Our results show that there is additional absorption in the 
red wing of Lyman-a that could be explained by accretion with velocities that are at least 
100 km s~^ and are likely as large as 290 km s~^. We envision a model in which a portion of 
the Lyman-a emission, likely formed near the accretion shocks, is reflected to the observer by 
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neutral hydrogen in the inner disk. A schematic cartoon of this model is shown in Figure ID 
The H2 pumping and fluorescence occurs where these Lyman-g photons ar e present in the 
inner disk, as described by the "thick disk" model in iHerczeg et al.l ( 12004 ). The reflected 
Lyman-a emission line, including the H2 absorption at the pumping wavelengths, is then 
absorbed by infalling neutral hydrogen in large accretion funnels in our line of sight. 

Accretion of gas from circur nstellar disks onto C TTSs is generally thought to be con- 
trolled by_stel]Bi;magiietic fields (IBouvier et al.l 120071 ). Strong magnetic fields of a few kilo- 
gauss (jJohns-KrullI 120071 . Yang et al. 2010, submitted) truncate the circumstellar gas disk 
at a few stellar radii and direct the accretion funnels onto the star near t he magnetic poles. 
From models of the spectral energy distribution (SED) of V4046 Sgr, iJensen &: Mathieu 
(I1997I ) found that dust in the inner regions is cleared out to about 0.18 AU, which is 38.6 
Rq. At this distance, the stellar magnetic fields are not strong enough to interact efficiently 
with the disk, and, more importantly, the temperature of molecular gas at the surface of the 
dusty disk may not be high enough for the hydrogen molecules to be electronically excited 
by the Lyman-a photons. For this transitional disk system, we think that both the fluo- 
rescent emission and accretion colurnns lik ely originate from the inner molecular gas disk 
(IMuzeroUe et al.l 120031 : iKastner et al.l 120081 ) that is closer to the central star than the dust 
disk. The absence of any significant differences between the radial velocities of the star, 
pumping lines and fiuorescent lines is consistent with the fiuorescent H2 gas lying in an inner 
gas disk seen nearly face-on. In this picture as shown in Figure HI w e expect to see more 
refiected light from the inner edge of the "thick disk" (cf. Figure 8 in lHerczeg et al.ll2004[ ) 
facing toward our line of sight than from the opposite side where the disk absorbs light along 
the line of sight to the observer. Therefore, the refiected Lyman-a emission line summed 
over all viewing angles toward both stars will show absorption mostly by downfiowing neu- 
tral hydrogen in the accretion columns (seen along the long black line in Figure H]) rather 
than upfiowing neutral hydrogen (seen along the long gray line in Figure H]) . 



Giinther fc SchmittI (120071 ) modeled the accretion shocks on V4046 Sgr, and their best-fit 
model of the X-ray observations yields an accretion rate of 3 x IO'^^Mq yr~^ and a maximum 
infall velocity of 535 km s~^ where the accretion gas strikes the stellar surface. Since the 
H I absorption velocity is between -1-100 and -1-290 km s~^, the a bsorbing H I gas is located 
between 2.3 and 4.4 R^, above the accretion shock (cf. Eq. 1 in lCalvet fc Gullbring||l998[ ) 
if we assume that the accretion columns have a constant cross-sectional area and the infalling 
gas sees only gravitational forces. 



This work was support by NASA grants NNX08AC146 and NAS5-98043 to the Univer- 
sity of Colorado at Boulder. The authors would like to thank the COS Science Team for 
help with data reduction and useful discussions. 
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Table 1. Lyman-a Absorption Features of V4046 Sgr. 



Wavelengths (A) 


Pumping Transitions 




Measured Eabs/Eem 


Optical Depth'' 


New Eabs/Eem 


1212.43 and 1214.54 


1-1 P(ll) and 1-1 R(12) 


0.0 and 0.0 


1.377 


0.055 


1.456 


1213.36 


3-1 P(14) 


0.015 


0.999 


0.087 


1.090 


1213.68 


4-2 R(12) 


0.050 


0.563 


0.106 


0.626 


1217.41 


4-0 P(19) 


0.417 


0.047 


3.505 


1.568 


1217.64 


0-2 R(l) 


0.0 


0.098 


1.700 


0.537 


1217.90 and 1219.98 


2-1 P(13) and 3-0 P(18) 


0.002 and 0.189 


0.226 


0.894 


0.553 



''Dissociation probability from upper level, calculated from lAbgrall et al. 1 l ll993h . 

''The optical depth is calculated for the additional absorption with Nhi = 2.5 X 10^^ cm~^ and v^^d = 290 km . 
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Fig. 1. — Two portions of the COS G130M spectra of DF Tau and V4046 Sgr. For display 
purposes, the DF Tau spectrum is scaled up by a factor of 4 and offset by +2 x 10"^'' erg s~^ 
cm~^ A~^. The red dashed vertical lines identify the fluorescent Lyman-band H2 emission 
lines and the red solid vertical lines identify the fluorescent Werner-band H2 lines. The 
strong emission feature around 1176 A is the C III multiplet. Two carbon monoxide (CO) 
emission features are also marked in the lower panel and will be presented by K. France et 
al. (2011, in preparation). 
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Fig. 2. — The Lyman-a emission profiles of DF Tau {top panel) and V4046 Sgr {bottom 
panel). The red dasfied lines indicate tfie uncertainty levels, and the vertical ticks mark the 
wavelengths of coincident Lyman-band H2 transitions. The blue areas indicate our estimates 
of the absorption in the H2 pumping transitions. In the bottom panel, the emission bump 
centered around 1215.7 A is geocoronal emission and is not from the star. 
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Fig. 3. — H2 absorption/emission ratios of DF Tau {top panel) and V4046 Sgr {bottom 
panel), corrected for additional Lyman-a absorption for a range of H I column densities and 
radial velocities. 
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Fig. 4. — A schematic cartoon of the V4046 Sgr system, which is not drawn to scale. 
The blue wavy arrows represent stellar and reflected Lyman-a photons, and the red wavy 
arrows repre sent the fluorescent H9. ph otons. The dust disk begins at 0.18 AU from the star, 
according to iJensen fc Mathieul (119971 ) . The short black arrows indicate the flow direction 
along the accretion funnels. The long black and gray arrows point to our line of sight, which 
is at ~35° with respect to the rotation axis of the disk. The long black arrow indicates that 
there is more reflected light from the inner edge of the gas disk that faces toward our line of 
sight than that from the opposite side, which is indicated by the long gray arrow, because 
of absorption through the disk along this line of sight. 



